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Abstract. 

We present a one-dimensional tight-binding chain of two-level systems coupled 
only through common dissipative Markovian reservoirs. This quantum chain 
can demonstrate anomalous thermodynamic behavior contradicting Fourier law. 
Population dynamics of individual systems of the chain is polynomial with the order 
determined by the initial state of the chain. The chain can simulate classically hard 
problems, such as multi-dimensional random walks. 
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1. Introduction 

Coherent chains of interacting quantum systems represent a wide class of physical 
objects that define the behavior of matter under different physical conditions. Study of 
theoretical models for these objects, which was started by the epoch-making works by 
Hubbard [I] and Ising [2], covers both new forms of matter and new types of interactions 
[3]. One can coherently chain cooled atoms in optical lattices [3], Josephson qubits 
in microwave transmission lines [1], semiconductor quantum dots, etc. Interactions 
between systems can be of quite different physical nature: spin-exchange interactions 
or pseudospin interactions corresponding to the dipole optical transitions [5], tunneling 
0, clipole-dipole interactions (3 [8j, photonic interactions (Jaynes-Cummings-Hubbard 
model m) and many others. Recent applications of these models led to a number of new 
fundamental results. For example, ID - chain of tunnel-coupled systems with two ends 
connected to heat reservoirs served as a model to justify the Fourier heat conduction law 
from the first principles and to provide microscopic definition of temperature pmn. 
Other results to be mentioned are the directivity of collective spontaneous emission 
HI EH), possibility to transfer quantum states [IB], the spatial propagation of Rabi 
oscillations (Rabi-waves) |CH, 15], and quantum optical nonreciprocity of the medium in 
timed Dicke state mt 

It is of particular importance to study new types of interactions that determine 
the coherent behavior of coupled systems. Here we suggest to couple them in a chain 
by connecting them pair-wisely to common dissipative reservoirs. It is already well- 
known that coupling several quantum systems to the same dissipative reservoir allows 
to obtain a number of highly non-trivial effects. First of all, such a coupling can create 
a decoherence-free subspace [18, 19]. An arbitrary initial state is eventually transferred 
into this subspace. Thus, it is possible to produce non-classical and entangled states 
via dissipative dynamics even without any direct interaction between quantum systems 
[20] . This effect has many potential applications. For example, it was recently used 
in practice to protect the highly-entangled initial polarization states of photons from 
dephasing in optical fibers PH- Moreover, combining coupling to the same reservoir with 
the nonlinear interaction between quantum systems, it is possible to produce nonlinear 
loss generating robustly a wide variety of non-classical states (221E3 EH- By adding 
external driving, a generated non-classical state can be preserved in the presence of an 
arbitrarily large linear loss [25]. 

Here we show that by coupling a set of quantum systems through common 
Markovian reservoirs and forming one-dimensional tight-binding chain, it is possible to 
produce highly non-trivial dynamics. Such a dissipatively coupled set of just a few two- 
level system renders a possibility to reproduce dynamics of far more complex systems, 
such as classical heat reservoirs and multidimensional random networks. This dynamics 
can be "tuned" by choice of the initial state of the system. 

We demonstrate that even for a few systems in the chain certain groups of 
matrix elements of the chain density matrix evolve according to equations formally 
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Figure 1 . (above) A scheme of bipartite tight-binding spin-chain with some elements 
subjected to strong loss. This scheme is described by Eq. (below) A scheme of 

dissipatively coupled tight-binding spin chain described by Eq. ([!]). 


coinciding with the equations describing classical random walks. For no more than 
one initial excitation in the chain, matrix elements of the single-excitation subspace 
evolve according to the equation describing two-dimensional classical random walk. 
Simultaneously, matrix elements describing coherences between the single-excitation and 
zero-excitation subspaces evolve according to the equation describing one-dimensional 
random walk. Taking more than one initial excitation in the chain, one is able to obtain 
dynamics governed by equations describing multi-dimensional classical random walks. 
This fact leads to a number of quite-counterintuitive consequences. Notwithstanding 
the Markovianity of reservoirs, a decay of excitation in any system of the chain is 
polynomial with the same power. By choosing the initial states, it is possible to obtain 
a population decay law l/t 2m+1 , for an arbitrary m > 0, provided that the chain is long 
enough. The sum of density matrix elements described by the random walk equation 
is conserved. So, one can observe rather peculiar thermodynamics-like behaviour of the 
chain. Despite seemingly classical character, this "thermodynamics" is quite anomalous. 
The stationary state of the chain can be entangled. Moreover, energy flow through the 
chain is not described by the Fourier law, whereas the flow of coherences is governed by 
it. 

The outline of paper is as follows. In the Section [2] we introduce the concept of 
tight-binding dissipatively coupled quantum chain. Also, we show how this "quantum 
gadget" can be built in practice from usual tight-binding unitary chain by selectively 
applying strong dissipation to certain systems in the chain. In the Section [3] we derive 
equations for the simplest case of the chain dynamics corresponding to an initial state 
with no more than just one excitation. Already in this case the chain exhibits a rich 
variety of highly unusual phenomena. We highlight the thermodynamical analogies 
in the Section |4| In the Section [5] we analyze in more details possibilities of having 
different polynomial decay for specific choices of initial states of the chain. Then, in the 
Section [6] we discuss the dynamics of the chain in the case of multiple initial excitations 
demonstrating how conservation of coherences might arise in this case, too. 
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Here we are discussing dynamics of a set of quantum objects coupled only through 
common dissipative reservoirs. We are interested in dynamics arising for the case when 
our set is rather large (with the number of objects, N 1). In this paper we consider 
one of the simplest cases of such a set, namely, a tight-binding one-dimensional chain 
of identical two-level systems (TLS), e. g., spins or pseudo-spins. The general scheme 
of such an arrangement is as follows. Each spin is supposed to be coupled to common 
Markovian reservoirs with one of two neighbors: jth TLS and (j + l)th TLS are coupled 
to common jth reservoir, (j — l)th TLS and jth TLS are coupled to (j — l)th reservoir, 
etc. Reservoirs are supposed to be independent (see Fig{l] (below)). Thus, we are 
considering here dynamics of the compound object described by the following generic 
effective master equation: 

d N 

Jt p = Y. v Kps; - ps+sj - sjsjp ). (i) 

ai j=1 

where Linbdlad operators are 

Sj — (Tj — aj +1 

and af = |±j) (=Fy I are rasing/lowering operators for the jth TLS. Vectors |±.,j 
describe excited/ground states of the jth TLS. Quantities 7 j are relaxation rates into 
corresponding reservoirs. For simplicity sake, we consider the finite-size homogeneous 
chain, thus we take 

7) 1 < j < N, 

0, j<0,j>TV+l. 

where 7 is a given constant and iV+1 is the total number of TLS in the chain, whereas N 
is the total number of reservoirs. The definition ([2]) corresponds to the total insulation 
of the chain from the outside areas of space, i.e., the absence of quantum flux over the 
chain boundaries. 

Notice that Eq. (JTj) is the quite general and describes a ID set of pairwise 
dissipatively coupled systems. A number of different physical models can lead to such an 
equation. For example, similar band structure of dissipators can arise in a graphenelike 
models based on a honeycomb lattice with nearest-neighbor and next nearest- neighbor 
couplings (Haldane model) [26]. Now let us discuss some simple cases as to demonstrate 
feasibility of dissipatively coupled chain and possibility to realize it in practice. 

First of all, let us point out that the dissipatively coupled chain can be produced just 
by subjecting some elements of a usual tight-binding chain with exchange interaction 
(like dipole-dipole or spin-spin interactions) to strong Markovian loss. Such a lattice with 
the regular pattern of dissipative sites (the every second one) were considered recently 
in the context of quantum walks and was shown to exhibit such non-trivial effects as 
topological transitions [33j. So, let us consider a simple example of a tight-binding 
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ID chain of TLS with every second TLS coupled to the separate bosonic reservoir (see 
Figjljabove)). Such a system is described by the following Hamiltonian 


Htotal (f ) Hq T Hchain “I” Hreseryoics T ^reservoirs i 


(3) 


where the Hamiltonian Hq describes non-interacting chain systems with the transition 
frequency Uq (we are using the system units with h = 1), 


H 0 = uj 0 J2 , ( 4 ) 

i 

and the Hamiltonian of direct spin-spin (or dipole-dipole) interaction is 

Hchain = W yX^IU-l a l. “ a f a l+l) ■ (3) 

l 


The Hamiltonian, H reservo i rs , describes modes of independent bosonic reservoirs coupled 
to the corresponding TLS of the chain 


H 


reservoirs 


= ^2^ikb] k bik, 


( 6 ) 


Ik 

where cis the frequency of the bosonic mode described by the annihilation operator, 
bik, and the creation operator b\ k . The operator V reservo i rs describes coupling of the chain 
systems to dissipative reservoirs, and can be represented as 


F reservoirs ^ 'X dlk&l bik + h.c.), (7) 

Ik 

where gik is the interaction constant for coupling of the Z-tli chain systems with the 
mode bi k . 

In the interaction picture with respect to the bosonic reservoirs Hamiltonian, 
Hreservoirs , and the chain Hamiltonian, H 0 , the total Hamiltonian (|3| becomes 

Htotal (t) = Hchain + V(t), (8) 

where 


V(t) = +h.c.), 


(9) 


and reservoir operators are 

= y vnbu, exp{-!(w,fc -wo)<}- 


( 10 ) 


We assume that our bosonic reservoirs are Markovian, mutually independent, and 
initially in the vacuum state. So, the following relations hold 

(R,(t)Rl(T)) K far ,s(t - t), 

(R,(t)R k (T)) = (R](t)Rl(T))=0, 


where T; is the decay rate into Z-th reservoir, 8(t — r) is the Dirac delta-function and 
the averaging denoted by (...) is carried over the states of reservoirs. 

Now let us suggest that only the every second system of the reservoir is subjected 
to loss, i.e. T 2m = 0. Also, we assume that losses are occurring on the time-scale much 
shorter than the dynamics of the excitation exchange described by the Hamiltonian (J5]) . 
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Then, taking for simplicity equal decay rates, T; = T, in the basis rotating with the 
frequency cuo, one arrives at the following master equation 

J t P = -i[Hchain , p] + E r ( 2 °2jP°H * P a tl a 2l ~ °tl°2lp) , l 11 ) 

Let us assume that TLS corresponding to dissipative sites are initially in the ground 
states. If the dissipation rate, T, is high enough in comparison with the strengths of 
direct interaction between neighboring spins, g, variables corresponding to TLS with 
even numbers can be adiabatically eliminated. Indeed, in this case for the dissipative 
sites we have 

= 0, ( cr 2l(t) a 2 m ( T )) = (^(^2m( r )) = °> 

(°' 2 ; 0 ) ( lL( t )) ~ hm exp{— 2 T(t - r)}. 

Then, taking TLS in dissipative sites as reservoirs and deriving the master equation 
up to the second order with respect to the ratio g 2 /T (see, for example, [33]), one can 
obtain from the master equation the master equation Eq. ([!]) describing dissipatively 
coupled tight-binding chain with the dissipative rate 7 ~ g 2 / 2 T. 

Thus, we have demonstrated one of the ways to obtain the dissipatively coupled 
chain. One can get it as the limiting case of a usual tight-binding chain with some 
sites subjected to losses. Such a bipartite dissipative lattice shown in FigjTJjb) can 
be realized in practice in a lot of different ways. For example, it can be build using 
a chain of color-center defects in diamond microcrystallites [35]. For color-center 
defects in diamond microcrystallites it is typical to have very low decoherence rates 
even at room temperature. Such objects can be manipulated with high precision and 
addressed individually by applied external electromagnetic fields [35]. Also, chains 
of individually manipulated atoms deposited on the metallic surface can be used for 
the purpose I37i m Eg. Among possible perspective candidates one can mentions 
schemes with trapped ions in optical lattices mm, photonic structures described by 
Jaynes-Cummings-Hubbard model [9], coupled system of optical waveguides pm. Bose- 
Einstein condensates in multiple-well potentials [22], or even networks of pigments in 
light-harvesting molecules mm- 

To realize the chain it is not necessary to manipulate "chain links" with high 
precision and address them individually. In our example, one can avoid coupling strongly 
dissipative reservoirs to the every second TLS in the chain. Restrictions imposed on 
precision of individual addressing can be significantly relaxed, if one considers a strongly 
coupled sufficiently long sub-chain instead of just one TLS as a lossy system to be 
adiabatically eliminated (an example of realistic consideration of such a dissipative sub¬ 
chain one can see, for example, in Ref. BZP- 

3. Single-excitation equations 

To clarify essential features of dynamics prescribed by the master equation ([!]), let us 
consider the initial chain state confined to the single-excitation subspace. Then, as 


Quantum tight-binding chains with dissipative coupling 


7 



Figure 2. (a) A solution of the master equation Eq. (121 for diagonal elements, pkk, 
for only one TLS initially excited (the one in the middle of the chain, j = 26). Dash- 
dotted, dotted and dashed curves correspond to the times yt = 3, 5, 9. Thin solid 


lines depict results Eq. (191 given by the heat equation, (b) A solution of Eq. (121 for 
all elements, pki, for only one TLS initially excited (j = 26); yt = 9. (c) Solutions 


of Eq.(12l for diagonal elements, pkk, for completely entangled initial states p± (241. 


Dash-dotted, dotted and dashed curves correspond to the times y t = 3,5,9 for the 
initial state p_, j = 25. Thin solid lines correspond to the same times for the initial 


state p+. For comparison, solutions are normalized, (d) A solution of Eq. (121 for all 
elements, pki, for the initial state j = 25; y t = 9. For all examples the number of 
TLS in the chain is N = 51. 


follows from Eq. 0, the single-excitation matrix elements, 

Pki = <l*|p|lz>, 

satisfy the following equation: 

d , . 

— — (7 k + lk-1 +7 Z + ll-l)Pkl + 

lkPk+ 1,1 + 7fc-lPfc-l,Z + 7zPfc,Z+l + 7Z-lPfc,Z-l5 (12) 

where the relaxation rates, 7 are defined in accordance with conditions (j 2 ]). Away from 


the edges, Eq. (12) coincides with a standard equation for a discrete classical random 


walk in two dimensions in continuous time, and describes diffusive propagation of the 
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excitation through the chain mediated by emission to reservoirs and re-absorption from 
them [28]. 

As follows from Eq. (JTj) , the coherences between the single-excitation subspace and 
the vacuum 


JV+l 


PkO — (lfc|p|0>, |0) — n | —A-) 


k =1 


satisfy the standard classical equation for the one-dimensional discrete classical random 
walk: 


rfjPkO — ~ (7*: + 1k-l)pk0 + lfkPk+ 1,0 + lfk-lPk-1,0' 


(13) 


The immediate consequence of the classical form of both Eqs. (12) and (13) is the 
conservation of sums of corresponding matrix elements: 

JV+l JV+l 

w = £ Pa (t) = y Pu ( 0 ), 

k,l =1 k,l =1 


7V+1 N+l 

F = ^2 Pko(t) = 5Z Pfco(0). 
k= 1 k= 1 


(14) 


Eqs. (14) are a consequence of collective character of coupling to reservoirs. The chain 


described by the master equation (jl| has another pure stationary state in addition to 
the usual lowest-energy one. This state satisfies the equation = 0, for all j. It is 

a pure and an entangled state. It describes the equal superposition of single-excitation 
states. For a chain with N + 1 elements this state is 


|T} = 


JV+l 


WTT S 


£ |lj>, 


(15) 


where 


ii?) - \+j) n 1 ~k)- 

k^3 


So, the state (15) describes just one excitation "spread" homogeneously over all the 
systems of the chain. 

Eqs. (12) and (13) can be easily solved analytically [33]. For example, the solution 


for the elements Pki(t) can be written as 


JV+l 


/>«(<)= £ “m» ex P {-+,+}. 


( 16 ) 


m,n=l 


where the eigenvalues for Eq. (12) are 


A m ,n = 4y |sin 2 + sin 2 


7T n 


N + l 


The coefficients at are defined from the initial state: 


a mn = VkmVln($m\p(Q)\®n), 
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where the coefficients are 


Vkm = \j( 2 - 5 k ,N+l) cos 


2nk 


-m 


N + l 

The symbol denotes Kronecker delta, and the vectors are 

1 


JV+1 


\$l) = 


'y ) ^inlln)- 


( 17 ) 


VN + 1 n=1 

Vectors 0 are mutually orthogonal, {$ TO |$ n ) = 5 mn . However, they do not diagonalize 
the density matrix, p. Since Ajv+i,jv+i = 0, the vector |$jv+i) is the stationary state 


|T), see Eq. (15). 


Solution of Eq.(13) is similar to the solution of Eq. (12). Eq. (12) can also be 


viewed as a discretization of the standard heat-transfer equation in the square L x L, 
where L = Na/2 is the chain length, and a is the doubled distance between neighboring 
systems in the original chain (the one depicted above in Fig. ([Tj)) . The later one can be 
obtained by assuming ak —)• x, al —> y: 


j t p{x ' V]t) 


a 2 y 


' d 2 p(x,y;t ) d 2 p(x, y, t) 


(18) 


dx 2 dy 2 J 

The absence of quantum flux over the chain boundaries corresponds to the Neumann 
boundary conditions for the density matrix at the square boundaries x = 0,L and 


y = 0, L. The discrete solution (16) approximate the solution for the Neumann boundary 
conditions only in the limit N —> oo. For initial states having only non-negative 
elements, pki(t = 0) > 0, Eq. (12) gives pki(t > 0) > 0. Thus, elements pki can 


be taken for classical probabilities, and our dissipatively coupled ID chain can indeed 
simulate 2D classical walk or heat transfer. In Fig.[2](a) examples of dynamics are shown 
for the case of just one TLS being initially completely excited. Even for modest number 
of TLS in the chain (N = 51 for our example depicted in Fig. (f2|) ) solutions of Eq. (12) 


given by Eq. (16) are very close to the fundamental solution of Eq. (18) obtained in the 


limit of inhnitely-length chain 
Po(x ’ v:t) = 

for 'jt > 1. Here Xo,yo denote the position of the initial excitation. 


exp 


{-[(* - x 0 ) 2 +(y- Sto) 2 ]/2a 2 7«} , 


(19) 


A continuous approximation for Eq.(13) is the ID heat-transfer equation, 


d . . 


2 d 2 p(z;t) 


( 20 ) 


All the consideration given above can be extended for this case, too. So, our dissipatively 
coupled ID chain can simulate simultaneously both ID and 2D classical random walk, 
or heat transfer. 

Non-exponential decay in our chain is stipulated by nonlocality of Lindblad 
operators in Eqs. (12) and (13)). It is interesting and instructive to compare the 


dynamics of population decay in our chain with dynamics of collective spontaneous 
emission of dense atomic cloud with dipole-dipole interaction mi- The continuous 
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limit of master equation correspondent to the model considered in Ref. HUE] is given 


by Eq.(A.5) in the Appendix A Eq.(A.5) looks similar to Eq. (18) and also manifests 


an algebraic law of population decay, 1 /y/jt [7]. However, it holds in 3D-space, and 
describes quite different physical process. The reason for it is the formal equivalence 


of Eq.(A.5) to the master equation with nonlocal Lindblad operators. This nonlocality 


arises from non-conservation of the excitation number and from accounting for quantum 
states corresponding to two excited atoms and one virtual photon with "negative" energy 

m- 

4 . Anomalous heat transfer 


Eqs. (12) and (13)) tell us that elements of the density matrix in the energy eigenstates 


basis evolve according to classical equations. However, despite this fact, the energy flow 
through the ID chain cannot be described by the ID classical Fourier law even in the 
limit of the large number of TLS in the chain. This holds even in the case, where the 
density matrix elements of the initial single-excitation state are non-negative. Indeed, 


as it was shown in the previous Section (see Eqs. (14)), the sum of matrix elements, 


N+l 

W = J2 Pkh 

k,l =1 

is preserved. Simultaneously, the total upper-state population of the chain does decay. 


As it follows from Eqs. (12), for example, for the initial state being just one completely 


excited TLS, the total population behaves in the following way: 

1 


I N+l \ N+l 

( H |lfc)(lfc| ) = 5Z Pkk 

\ k =1 / k=1 


OC 


Vrt' 


A physical reason for breaking the classical Fourier heat conductivity law is rather 
simple. Besides the energy transport occurring in one dimension (along the chain), 
there is an additional motion. Fourier law is always breaking down in presence of 
additional motions (the simplest example is the convection in liquids |29|). In our 
case this additional motion is a flow of quantum coherence. Indeed, in the continuous 
approximation the energy flux along the chain reads as 


d 

J(x,x ) oc lim — p(x, y, t). 
y ’ ’ v^dx y ’ 


( 21 ) 


The right-hand part of this equation cannot be represented in the form of the gradient 
of a scalar function satisfying the ID heat-transfer equation. Thus, the right-hand 


part of Eq.(21), or the energy per TLS, p(x,x;t), cannot be associated with classical 
Since in our chain not only energy is transported, but also quantum 


temperature 

coherences, it is possible to introduce the concept of the effective 2D "quantum 
temperature" and 2D quantum flux as, correspondingly, 

T q (x, y\ t) = p(x,y\t ), 

Jq(x,y,t ) oc V xy p(x,y,t). 
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The value of the real energy flux corresponds to the component of quantum flux, 
e x • J q (x,x‘,t ), where e x is the unit vector along the chain, while the ^/-component 
describes the coherences transfer. This additional motion produces a second dimension 
in the heat transfer equation and leads to the anomalous thermodynamics. A number 
of examples of an anomalous thermodynamical behavior for quantum structures is 
known (see, for example, Refs. OH EDGE]). But the dissipatively coupled chain differs 
starkly by the nature of anomaly. One can retrieve classical thermodynamics and 
establish connection between "classical" and "quantum" temperatures by neglecting 
the additional motions in Eq.(18), i.e. by averaging out some coherences. Indeed, 


introducing the averaged temperature and averaged energy flux as 


T, 


class 


(x;f) = J p(x,y\t)dy, 


Jcim(x\t)= / J q (x,y;t)dy 


( 22 ) 

(23) 


one obtains ID heat transfer equation for temperature and relation between the 
temperature and energy flux, 

-» d 

^class(*D bx ^ -^"class(*D f: 

in full agreement with Fourier heat conductivity theory. But there is no convincing 


reasons for such approximation to be done and for the association of the integrals (22 


23) with a real physical temperature and a real energy flux. As it will be shown below, 
use of such classical model in our case leads to the disappearance of some important 
observable regimes. 


5. Polynomial decay of populations 


The dissipatively coupled chain offers unique possibility of controlling the quantum 
state dynamics. Just by choosing different initial state, one can manipulate the law of 
spontaneous decay. Also, dependence of the decay law on the type of initial state gives 
one an opportunity of identifying initial states of the chain by measuring the population 
of just one TLS. 

We have obtained that for the simple ID set of TLS coupled to common Markovian 
reservoirs and the single initially excited TLS, the upper-state population decay for 
chain TLS occurs polynomially, according the the 1/t law. Let us now demonstrate 
by just adjusting the initial state of the chain, one can obtain a population decay law 
1/t 2,71+1 for an arbitrary m > 0. Interestingly, we can arrange initial conditions that are 
hardly possible in the classical case. In particular, it is possible to create equivalents of 
neighboring regions with "positive" and "negative" temperatures (or even "imaginary" 
ones). It can be done by choosing an initial state with components orthogonal to the 


stationary state (15). An example of pu dynamics for such a state is demonstrated in 


Fig. He ), where the solutions of Eq. ([I2| are shown for entangled initial states: 


P± = IV’iXV’: 


±1 


\^±) = -^(i 1 ?) ± iw))- 


(24) 
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Figure 3. Contour plots of exact solutions of the master equation Eq. (12) for 
initial states ( 25 1. For the panel (a) the initial state is the just one completely excited 
system, m = 0. For the panel (b) m = 1, /o = fi = 1/V2. For the panel (c) 
m = 2, / 0 = f 2 = l/x/6, /i = 2/\/6. For the panel (d) m = 3, f 0 = f 3 = l/v^O, 
— 3/a/ 20- For all the panels 'yt = 9, N = 51. For comparison, solutions are 
normalized (the largest element is taken to be unity). 


The solution for the state p + behaves itself "classically": both closely situated initial 
peaks are soon merged together in one Gaussian shape. However, dynamics with the 
initial state is drastically different: negative elements, pki, do persist (Fig. 1(d)), 
and initially close peaks are distancing from each other with time. Moreover, excitation 
displacement for the initial state p_ appears to be faster than for the initial state p + . 
Such a difference can be quite simply illustrated with the solution of the continuous 
approximation (18). The continuous analog of p± can be represented as 

P±(x, y\t = 0 ) oc 5(x - x 0 )5(y - y 0 ) + S(x - (x 0 + a))5(y - (y 0 + a)) 

±5(x ~ (x 0 + a))5(y - y 0 ) ± 5(x - x 0 )5(y - (y 0 + a)), 

where 5(x) is the delta-function. It is easy to see that when the excitation has spread 
for distances x, y such that 

\j{x- x 0 ) 2 + (y- y 0 ) 2 > a , 


the solution for the initial state p + is given by Eq. (19). However, for the initial state 
P- dynamics is quite different. When the excitation has spread far enough, the solution 
for the initial state p_ is 

p{x,y\t) « 2p 0 (x,y;t)(x - x 0 )(y - y 0 )/(2a^yt) 3 . 

Thus, just by choosing the initial state to be p_, we get the population decay law 1/f 3 . 
This state doesn’t exist within the bounds of quasiclassical Fourier theory. Indeed, the 
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averaging of this state following (22, 23) leads to the trivial solution: the temperature 


as well as energy flux equal to zero at arbitrary point of space for every moment of time. 
It is not hard to see that the population decay law can be further varied by choice of 
the initial state. Let us consider pure initial states of the general form 


'-k+i 


>, 


(25) 


W = B-t'/illi 

i =0 

where fi are some scalar coefficients, and the excitation is assumed to be well localized, 
m <C N. Taking the continuous approximation and using the quantity m/'yt C 1 as 


small parameter, from Eq. (25) we straightforwardly obtain the following approximation: 

1 


p(x,y;t) « po(x, y, t) £ 


n = 0 n\(2a 2 jty 


*£( 

hj =0 


- 1 ) 


i+j 


fifj[i( x - oco) + j{y 


y 0 )Y 


(26) 


Starting with the initial state (25) and choosing the coefficients fi in the form of the 


normalized binomial coefficients, fi oc m\/i\(m — i)\, we get the first non-zero term in 


the approximation (26) of the order of l/t 2m+1 . Examples of exact solutions for initial 


states (25) are shown in Figs. [3] The panels (a-d) corresponds to m — 0,1, 2, 3. 

Notice, that the states (|25| are non-correctly described within the bounds of the 


quasiclassical Fourier theory. The "classical" averaging of states (25) as it is done in 


Eqs. (22, 23) can lead to the drastic deformation of them (for example, the states 


orthogonal to the stationary state will be reduced to zero by such averaging). Mind 
that polynomial decay regimes considered here are taking place for times, when initially 
localized excitation has spread far enough. However, we do not consider here an influence 
of edges. That is, we limit ourselves to the time intervals satisfying am < v(t) < aN, 
where v(t) is the variance of the distribution pki(t). For example, for the single initially 
excited TLS this condition reads as 1 < yf < N. An influence of chain edges will be 
considered in further work. 


6. Multi-excitation dynamics 

It is easy to surmise from Eq. ([Tj) that there is a profound difference between chain 
dynamics in cases of one and multiple initial excitations. Indeed, the state with 


multiple excitations will inevitable decay toward mixture of the stationary state (15) 


and the vacuum state. Generally, for the initial state with no more than N excitation, 
equations of density elements corresponding to the IV-excitation subspace describe 2N- 
dimensional classical random walk with losses (see equations for an arbitrary number of 


initial excitations given in Appendix B). However, even in this case one can still model 


lossless multi-dimensional random walks within some limited period of time and for 
some specific classes of initial states. Moreover, there are regimes when one can have 
preservation of coherences, where the chain behaves itself in quite "thermodynamic" 
way. 
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Let us illustrate our consideration with the case of no more than two initial 
excitations in the chain . Two-excitation matrix elements are 


where 




if) — +fc)i+i) n i —})■ 

j^k,i 


In difference with Eqs.(12), there are two different kinds of equations for the matrix 
elements p^n- For the elements without neighbouring indexes, i.e. for k 7 ^ l ± 1 and 
m 7 ^ n ± 1, from Eq.Q one gets the following system of equations (for the sake of 
illustration we are giving here equations only for the internal TLS of the chain, i.e. 
k, l, m , n 7 ^ 1, IV + 1): 

1 d 


7 dt 


p k ’ 1 — 
rm,n 


= —8 n k ’ 1 4- n k+1 ' 1 -I- n k ~ 1,1 
°rm.n ' rmji ' rn 


'm,n 


I fc,/—1 I k,l+l I k,l 
' rm,n ' rm,n 1 rm—l,n 
■ k,l - k,l - k,l 
' Pm+l,n "t” Pm,n— 1 ' Pm,n-\- 1* 


(27) 


Obviously, Eq.(27) coincides with the equation for 4D classical random walk and in the 


continuous limit transforms to the 4D heat-transfer equation. 

The situation is quite different for the matrix elements with neighbouring indexes. 


Let us assume, for example, that l — k + 1. Then, instead of Eqs. (27) one has 

1 __ o k,k+l , k—l,k+l , k,k +2 

L Pm,n ®Pm,n ' Pm,n ' Pm,n 


7 dt 


■ k,l . k,l 

'Pm—l,n ' Pm+1; 


■ k,l - k,l 

' Pm,n—1 ' Pm,n+1* 


(28) 


Eq.(28) does not coincide with the equation for classical random walk. It contains terms 
describing loss. It can be easily seen in the continuous limit, since then, for example, 
one has 


2 

Pm.n 


k,k+i 
rm.n 


.. d . . 

hm a-—p(xi,x 2 ,y,z;t), 


X2—yzi Qx2 

where variables X\,x 2 ,y,z correspond to the indexes k,l,m,n, respectively. Of course, 
it should be expected. Multi-excitation state does decay toward the single-excitation 
one. 


Eqs. (27) and (28) point to a number of quite counter-intuitive conclusions. First of 


all, if the initially excited TLS are far from each other, the spread of coherences occurs 
as for lossless multi-dimensional random walk till the excitation spreads to neighbouring 
TLS. Notice, that it will take place even in the case of entangled initial state of several 
TLS. If one has initially several uncorrelated excited TLS being far from another, the 
whole chain behaves like a set of unconnected chains with just one excitation per chain 
till the excitation spreads to neighbouring TLS. So, it means that one might have a sort 
of conservation of coherences in this case, too. 

Let us demonstrate an appearance of such a "conservation law of a sorts" with 
the example of two initially non-interacting chains with no more that one excitation in 
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each. So, we take that there are chains of M and N TLS with both M,N 1. Also, 
corresponding sums of coherences (14), denoted as Wm and Wx are supposed to be 
much less than number of TLS, Wx X, X = M, N. We assume that the chains are 
initially in the stationary states, 


/# = 


W x 

X 


+ 1 


Wx 

X 


\V X ){V> 


A'| 


(29) 


where X = M,N , the vectors |\Ew) are given by Eq.(15) and vectors \Vx) describe the 
lowest energy state of corresponding chains. 

Now let us assume that the chains are coupled by the common dissipative reservoir 
connecting the last TLS of the first chain and the first TLS of the second chain (the rate 
of decay to this reservoir we take to be the same 7 as the rates for all other reservoirs). 
Then, the stationary state of compound chain will be obviously given also by Eq.(29) 
with X = M + N. ft is easy to get from Eqs.(27) and (28) that the sum of coherences 
for the new stationary state is 


Wn+m — Wx + Wm + O (^WxW m A yy)) ~ Wx + Wm- (30) 

Notice that the expression ( |30| becomes exact in the limit M, N —» 00 . So, the sum of 
coherences for the stationary state of the compound chain is indeed approximately equal 
to the same of coherences of parts. Moreover, if one disconnects chains and they settle 
into stationary states again, the value of matrix elements of single-excitation density 
matrix (the "quantum temperature" as introduced in the Section [4]) for both part will 
remain to be equal, = (Wm+n/(M + N ) 2 ). 


7. Conclusions 

We have suggested and discussed a tight-binding dissipatively coupled quantum chain. 
It consists of two-level systems pairwise coupled to the same Markovian reservoirs. We 
have shown that despite being composed of a comparatively few systems forming the 
simplest one-dimensional chain, such a chain can model a multi-dimensional random 
walk, or model a solution of multi-dimensional heat-transfer equation. In difference 
to the classical "original", it is possible to model heat flow from initial distributions 
with regions of positive and negative temperatures. Such possibility of quantum 
stimulations in many-body physics corresponds to the long-standing idea first proposed 
by R. Feinman and recently cited by C. Cohen-Tannoudji and D. Guery-Odelin as a 
conclusive remark to their book (3]. This model exhibits anomalous thermodynamics 
behavior, such as non-Fourier heat conductivity and non-existence of temperature in 
the classical meaning. The population dynamics of all TLS in the chain is always non¬ 
exponential and exhibit polynomial character. By the choice of the initial state, different 
power laws of decay, l/t 2m+1 , for an arbitrary m > 0 can be achieved. The suggested 
chain can be used as an efficient simulator of classically hard problems, such as multi¬ 
dimensional quantum walks or heat-transfer equations. Note, that the considered chain 
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is not unique. One can build a number of different dissipatively coupled chains, for 
example, by changing phases in the system-system interaction terms in the chain. 
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Appendix A. Comparison of dissipatevely coupled quantum chain with 
spontaneously emitted atomic cloud 


It is interesting and constructive to compare single-excitation dynamics of the 
dissipatively coupled quantum chain with collective spontaneous emission in the dense 
atomic cloud with dipole-dipole interactions (see Ref. 0 ). In 0 in the main text a system 
of N TLS is considered. Initially, one of them is in the excited state, while all others 
are in the ground state. TLS are placed at positions r}, j = 1,2... A. Multi-mode 
electromagnetic held is interacting with all TLS. Initially, this held is in the vacuum 
state, | vac). The solution of the Schrodinger equation for the atoms and held can be 
represented as 


N 

W)> ='52Pj( t )\ l i)\ VaC ) + \^rest(t)), (A.l) 

3= 1 

where the coefficients 8j(t) are amplitudes of having jth TLS completely excited; the 
wave function 'I' rest (f)) denote components of the total wave function, |T(f)), having 
other then single excited TLS and the held vacuum. Notice that in Ref. [ 8 | of the main 
text the rotating-wave approximation was not used to describe TLS-held interaction. So, 
the function |^/ rest (f)) includes also components with the number of excitations larger 
than one. 

Using standard approximations about the reservoir, in Ref. [ 8 ] the following system 
of equations was obtained 


4&(f) = -7 Pj(t) + ?7 Pk(t)K jk , 

k^j 


Kjk ~ 


exp{-ik 0 \r k - fj\} 
k 0 \r k - fj | 


(A. 2 ) 


(A.3) 


where k 0 = ooq/c, ujq is the TLS transition frequency, 7 is the decay rate into the held 
reservoir. 


Eq.(A.2) can be easily transformed to the form similar to Eq.(12) of the main text. 
Indeed, introducing the matrix elements pj k = f3j/3 k , from Eq. (A.2) it follows that 

xrPjfc = Pj k + 27 ( PiiKji ~ Pik^ik ) • 


Continuous analog of Eq. (A.4) reads 

r) = ~2lp{r, r) 


+ i ~i^j< eR 

v 

-hy ji 3 R 

V 


exp{— ik 0 \r — i?|} 
k 0 \r — R| 

exp{f/co|r — R\} 
k 0 \r — R\ 


P(R, r) 


p(r,R), 


(A.4) 


where V is the volume of the cloud. Equations Eq.(A.4), Eq.(A.5) are similar to the Eqs. 
(12), (26), respectively. The difference consists in the form of operators in right-hand 
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parts (integral operators with the Green function of Helmholz equation in the capacity of 


kernel in Eq.( A.5) instead of Laplace operator in Eq.(26). Integral operators in Eq.( A.5) 


are non-Hermitian due to the presence of counter-rotating terms in Eq.(A.2) which are 


stipulated by the virtual photons with "negative" energy. It lead to the matrix elements 


(A.3) instead of 


hj k i 


.sin{7c 0 |f fc -fj |} 


ko\r k - fj 


(A.5) 


As a result, Eq (A.4) is not of Liouville type. It is equal (as Eq. (26) in the main text) 


to the master equation with non-local Linblad operators. However, it describes rather 
different physical process. Thus, as it was shown in Ref.[H] in the main text, collective 
spontaneous emission manifests for some types of boundary conditions algebraic law of 
population decay. Such behavior is stipulated by the virtual photons with "negative" 


energy and disappears in the rotation-wave approximation (done by replacing Eq. (A.3) 


with Eq.(A.5)). 


Appendix B. Equations for the matrix elements for the case of multiple 
excitations 

Here we write down equations for matrix elements corresponding to states with the 
highest possible number of excitations. Let us assume that initially we have no more 
than m excitations in the chain (for example, m TLS in completely excited state). From 
Eq.Q it is not hard to obtain the following set of equations for the matrix elements 
corresponding to m-excitation subspace 


d N+l 

~r:{ m ki,k2,-.,k ni \p\ m h,i2,...,im.) = A Ah) 

ai i=1 


(B.l) 


where 


~ (m {k . } \)}p\m {li} ), (B.2) 

Pi = {m{k i }\p{'Yi i (\ m {i j +6j,i}) ~ l m (hl)) 

+7ii-i(l m {k-Ad> _ I 

Here the state |mp.}) = is the state with m TLS with numbers 

being completely excited. Also, \m{i j±Sj4 y) = As it is 

in Eq.(3) of the main body of the paper, we are assuming that y 0 = 0, and 7 fc >0 = 7 . 
Notice that in case when the indexes are neighbors, i.e. when k i+l — A ; 8 + 1, one has 

Ai = —7fei( m {fej}|p| m {ii}) + 

7 fc i -i(^{fc J -5 i , i }|p|m {Zi} ) - (m {kj }\p\m {li} ), 

Ai+i = 7fej+i ( m {kj+5j, i+ i}\p\ m {ii}) ~ 

lk i+1 (m {kj} \p\m { i i} ) - 7 k i ( m { k j} \p\m { i i} ). 
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The same holds for the case of neighboring h and li+\. For indexes equal to 1 or N + 1 
corresponding As are 

“ 7fci((' m fci+l,fc 2 ,...,fc m | — ( m k 1 ,k 2 ,...,km\)\p\ m {li}) i 

\ N +1 = - ( m kl,k2,-,km\)\p\ rn {h})- 

For matrix elements without neighboring excitations, i.e. ki ^ kj ± 1, U ^ lj ± 1, 
\/ki,kj, one obtains that Eq.(B.l) is the discretization of the following 2m-dimensional 
heat-transfer equation 


d (f xs x 2 ^ d 2 p({x k },{y k }-,t) 

w, 1 Vkh o~ a 7L- 


dx\ 


+« 2 7 E 


k =1 


d 2 p({x k },{y k }-,t) 

dyl 


(B.3) 






